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O' 

O ' Abstract 

(n: 

5_^ . In relativistic heavy ion collisions an exact multipole decomposition of the 

Qh| Lorentz transformed time dependent Coulomb potentials in a coordinate system 

■^ ' with equal constant, but opposite velocities of the ions, is obtained for both zero 

cn ' and different from zero impact parameter. The case of large values of 7 and the 

'"^ , gauge transformation of the interaction removing both the 7 dependence and 

,__! ' the In h dependence are also considered. 
>■ 

§ ; 1 Introduction 

^: 

\^ • Relativistic heavy ion collisions offer us the opportunity to investigate the properties 

^ \ of nuclear matter under very extreme conditions of energy and temperature [1]. In- 

^— > i formation on the collision processes is only obtained indirectly from the energy and 

Oh. angular distributions of the produced particles and from fragment distributions. Elec- 

,-S^ ' tromagnetic probes like the photon and dilepton can yield esential information on 

. ^ \ the evolution of relativistic heavy ion collisions, since these signals leave the reaction 

rS \ zone without distorsion by strong interaction. Photon and dilepton production are 

traditionally considered to be one of the possible probes for detecting the quark-gluon 
plasma [2]. 

Relativistic heavy ion collisions provide a tool for the investigation of electrons in 
extremely strong electromagnetic fields [3]. Different methods were used for studying 
the atomic physics effects - coupled channel equations, perturbation theory [4], finite 
difference equations (electron-positron pair production, pair production with capture). 
The field of electron-positron pair creation in relativistic heavy ion collisions is reviwed 
in Refs. [4, 5]. Different processes were studied in Refs. [6-9]. Electron-positron pair 
creation in relativistic heavy ion collisions was non-perturbatively described by coupled 
channel equations in a coordinate system with equal constant, but opposite velocities 
of the ions [10]. Electron-positron pairs are produced by the electromagnetic fields 
of high Z-ions in atomic collisions at relativistic velocities. Atomic collisions mean 
collisions with impact parameters which are larger than the sum of the nuclear radii. 



c^ 



i.e., the nuclear forces have no influence on the collision. The problem of dileptons 
excited through the mechanism of virtual bremsstrahlung was investigated in Ref. [11]. 
In Sec. II we examine the basic step of the calculational technology: the multipole 
expansion. An expansion in I/7 produces the multipole expansion as compact, explicit 
and simple functions of space and time. We will fix the coordinate system in the 
centrum of mass of the colliding ions, which is the more appropriate system for studying 
the multiple pair production process. We proceed to a multipole decomposition of the 
two potentials, in order to approximate the two-center Dirac equation by a simpler one- 
center Dirac equation. We shall consider both the cases of large values of 7 and exact 
multipole decomposition, for both different from zero and zero impact parameter. In 
Sec. Ill we discuss about the gauge transformation of the interaction removing both 
the 7 dependence and the In 6 dependence. 

2 Classical form of the interaction 

In the consideration of the electromagnetic interaction of relativistic heavy ions, their 
motion is assumed to be well described by straight line, classical paths and unperturbed 
by recoil effects. Since we are interested in the multiple pair production process, it 
is more natural to fix the coordinate system in the centrum of mass (or center of 
momentum - we consider the same type of colliding ions, with the same A and Z 
and equal constant, but opposite velocities) of the colliding ions, rather than on one 
of the ions, which is the case, for example, for bound electron-positron creation. The 
interaction of the electron-positron field with the moving ions is then given by the 
Lorentz transformed Coulomb potentials (Lienard-Wiechert potentials) 

... ,. eZ(l-va) 

Vi{p,z,t) 



{[(b/2-p)/7]^ + (z-^t)2}V2' 

gZ{1 + vaz) 
V,[p, z, t) = ^^^_^^2 _ p)/^]2 + (^ + ^tyy/2 ■ (2) 

Here b is the distance between the ion straight line paths, which are taken along the z 
axis (see Fig. 1); eZ are the charges and v, — v the velocities of the ions, 7 = 1/(1 — 
f ^)^/^, p, z and t are the coordinates of the electron-positron field relative to the cm. 
system and a^ is the Dirac matrix. The characteristic effects of the electromagnetic 
interactions eVi and eV2 produced by the ions are contained in their time dependence 
and the very severe compression of the spatial dependences (equivalently - the sharp 



pulse description). The time-dependent two-center Dirac equation for the electron- 
positron motion is given by 

i^^^{r,t) = HD^{r,t) (3) 

(we employ the natural unit system with h = c = rrie = 1). The Dirac Hamiltonian 
Hd is given by 

Hd = Ho- eVi - eV2, (4) 

Ho = c-p + p. (5) 

Since a method for constructing continuum wave functions for the two-center 
Dirac equation is known only for the static case (only for scalar potentials) [12, 13], 
we proceed to a multipole decomposition of Vi and V2, in order to approximate the 
two-center Dirac equation (3) by a simpler one-center Dirac equation. We will consider 
the following cases: 

A. - large values of 7 and 6 7^ (but not too large); 

B. - exact multipole decomposition and 6 7^ 0; 
C - exact multipole decomposition and 6 = 0; 
D. - large values of 7 and 6 = 0. 

2.1 A. The case of large values of 7 and 6 7^ (but not too 
large) 

In [14] there were derived simple closed forms for the multipole M™ that are accurate 
for large values of 7 (asymptotic-large 7 limit) for the transformed Coulomb potential 
Vp in the case of a moving projectile and the coordinate system fixed on the target 
nucleus: 

V{pzt) = e^p(l - Vp^.) .g^ 

where Zp, v^ and b are respectively the charge, the velocity and the impact parameter 
of the projectile. This multipole decomposition of the potential Vy reads 

unr, t) ^ J ''"^• {|(b-P)/of + (.-v)'}'/- <" 

VM z, t) = eZ,{l - v,a,) Y, Mr{r, t)Yr{l>, <l>) = eZ,{l - v,a,) Y, V,"'(r, *)■ (8) 

Lm Lm 



The derivation consists of an almost straightforward expansion about the near singular 
point z = Vpt, after some careful rearrangements of the integrand. It is important to 
note that the expansion requires that 6/(r — Vpt) ^ 7; assuming that the nature of the 
physical problem confines r and Vpt to moderate multiples of h/rrieC, the requirement 
of b becomes b ^ 'yh/rrieC [14]. Defining Of by 



(1) \vpt\ < r, 



cos I 



^-, sin 6^4 
r 



^i-Cff, 



(2) \vpt\ > r, cosOt = 1, sin 6*4 = 0, 



(9) 



replacing Vp by 1 (they differ only by 0(1/7^) and noting that to order In 7/7^ the 
(b — p)^/7^ term in the denominator can be replaced by its value at ^ = 6t, the 
asymptotic form of the multipole decomposition of Vp{p, z, t) (where z = r cos 6) is 
given for m > by [14, 15] 



Mr{r,t) 



yrii^) 


27r r2-t2 


r 


m^ bl 
27r b^ 



0, 



r <t 



t <r < VP + t^ 



(10) 



> V^Tt^. 



This expression is valid for positive m and t. Negative m and t expressions are given 
by symmetry. For m = 0, the asymptotic form is [14, 15] 



M^{r,t) 



\/^(2/ + 1: 



X < 



2Qi 



r <t 



Pz(-)(21n27 + ln^ 



t 



e 



-25:-), t<r<v/FTI^(ii) 



n=l 
I 



PK-)(21n27-2^ 



n=l 



>.JW+¥, 



valid for positive t, with the negative t form again to be obtained by symmetry. Here 
Qi is the Legendre function of the second kind. Unless the m = multipole terms are 
to be integrated together with functions that vary sharply in the r ~ Vpt region, the 
need for the bridge function at r ~ t disappears [14]. The evaluation of time-dependent 
matrix elements for use in coupled-channels calculations will be facilitated by the form 
of the above interaction (which will have an analogous structure for our two-center 
case). In Eqs. (10), (11) all terms can be expressed in terms of polynomials and/or 
convergent series in t/r for t < r and in r/t for r <t. Thus the matrix element to be 
integrated over r can have its interaction expressed as a series of negative and positive 
powers of r. In each of the terms of the series in r, t is effectivelly a coefficient: r and t 



have become separable term by term. These asymptotic forms provide a simple means 
of calculating matrix elements. In addition they allow to straightforwardly look at the 
7 dependence of any set of matrix elements or amplitudes. As one can see directly, the 
7 dependence of the multipole operators appears only in the m = terms and only in 
the one form [14] 



2{\n2^)eZ,{l - a,)Y^Yi\9) 



7r{2l 



r 



2(ln27)eZp(l - a,)5{z -t), r> t;(12) 



the condition r > t is automatically contained in the 6{z — t). But as we shall see in 
Sec. Ill, this set of terms is entirely removable by a gauge transformation. Once having 
removed this explicit dependence on 7, we are left only with the implicit dependence 
inherent in the limitation already noted that the above forms are valid only for the 
region b <^ '-fh/meC. However, when b is not small, perturbation theory is valid. 

The multipole expansion for the potential Vi is given also by the formulas (10), (11), 
but with b replaced by 6/2, Zp by Z and Vp by v with v = 1. We have proven that the 
multipole expansion for V2 can be obtained from that for Vi, by replacing 1 — a^ by 
1 + ^2 and by multiplying the whole expression by (—1)'. 

1. Monopole approximation 

In the case / = 0, m = 0, we have 



w 



(0) 



2Qc 



.t. 



r <t 



eZil — Uz) — X < 



21n27 + ln 



6V4 ' 
2 In 27, 



t <r < 

r > Jb^/4: 



^b^/A + 12 (13) 



t2 



-(0) _ T/{0) 



and Vr = Vr{a, -* -a,). Then 



V(°) = Vi'^ + yf = ^ X < 



Qo(-), r<t 

r 

ln27 + -ln^^^^, t<r< ^b^A + t^ 



(14) 



In 27, r>J674 + t2. 



Here 



r 2 t — r 
2. Dipole approximation (I = 1) 



(15) 



The dipole term is given by 



\/(i) = vl'^ + V,^'^ = -^^^ 



r <t 



[sin 6 cos (J) sin 6t x < 



' by 4 



t <r < Jb'^/4: + t^ 



.2_^2-)^ r>y'6V4 + t2 



Here 



+ cos 61 X < 



2gi(-), r<t 

r 

^(21n27 + ln^— ^-2), t < r < ,J¥j4 + ¥ ]. 



(16) 



2-(ln27-l), r>y'6V4 + t2. 



Qi(-) = -go(-)-l = f ln^-1. 
r J- J- 2r t — r 



;i7) 



We remind that f = 1 in the case A. In the dipole approximation we have to add the 
terms V^^^ and V"'-"'^^ given by Eqs. (14) and (16): 

\/ = \/(o) + \/(i). (18) 



2.2 B. Exact multipole decomposition (6 7^ 0) 

The exact multipole decomposition is given as the series [14]: 



1, 



Vi = eZ(l - va,) J2 YnO, 0) J2 y^^iK 0)7^(/, /'; r, u)A{m- /, /'; -), (19) 



where 



and 



n{l,l';r,u) 



r / \2il /2 n ^ • n 

u=[— + {vt)\^, cos9u = —, sm6'„ = -- 
4 M 2m 



TT W 



r 



■;+;'+!' 



(20) 



4r''+ir(^ + i)r(/' + 



^J + /' + l /'-/ ,, 3 m2 

-i^f = ^—J+7^,Z^) for r>M,(21) 



2 ' 2 



2'r^ 



A{m;l,l';^) = ^(-l)™v^(2/ + 1)(2/' + 1) EC"'^^oC^1^!Ql(^- 



(22) 



For r < -u, we have to do the interchange r ^^ u, I ^^ I' in TZ. For the special case 
r = M, 



vrl 



r 



■i+i'+i^ 



^^^' ^'' ""' ""^ 4 r r(^ + i)r(^ + i)r(^ + 1) ' 

6 



(23) 



This form for the muhipole expansion is equivalent with that appearing in the hter- 
ature. There is a set of special cases that is particularly simple [14]. For r > u the 
sum over /' in Eq. (19) is cut off by the spherical Bessel function integral (for TZ) (the 
occurence of a negative-integer value of the argument of the gamma function ensures 
that I' < I); furthermore, for /' < / the hypergeometric function is just a polynomial. 
Thus the exact result for Y"/ component of Vi is for r > u: 



Y,\d,, 



/Stt b 
2 2^ 



V V V V 



(24) 



^Qo(^- 



(25) 



As 7 — s> CX3, w — !► 1, then the term in square brackets approaches —1, in agreement with 
the corresponding expression in the case A. Similarly, the Y^ component of Vi is for 
r > u: 

r Iv V 

As 7 -^ oo, the term in square brackets approaches 2(ln27 — 1), in agreement with 
the case A. The difference between exact and asymptotic values is O (In 7/7^). 

Another form of the exact multipole decomposition can be obtained if we write 
the denominator of Vi as 

1 1 



^{x-b/2y + y^ + -f^{z-vty 



R'(t)|' 



(26) 



with the vectors r' = {x,y,'~fz) and R' = (6/2, 0,7ft). The expansion in multipoles is 
then given by [3, 16]: 



R' 



°° 47r 

E 
1=0 

00 

E 



2/ + 1 R'^ 
An R 



- y: Yr{R')Yi 



m/^'\ 



m=—l 
// I 



1=0 



21 + 1 r"+i 



E yr{R')yi 



""r'] 



r' <R' 
r' > R'. 



(27) 



m=~l 



R' and f' represent the angular arguments of the spherical harmonics Y"^ which give 
the direction of R' and r'. As usual the scattering plane is taken to be = 0. In 
consequence the spherical harmonics Yi"^*{R') is real- valued. 
1. Monopole approximation {I = 0) 



We have (with 7^^^ =7^ — 1 and R' = Jb'^/4: + 'y'^vH'^/'y): 



w 



(0) 



eZ(l — va^) X 




(28) 



and 14^°^ = Vl°\a, -^ -a,). Then 



, — , r'<R' 
\/(°)=2eZx H' ~ (29) 



2. Dipole approximation 
For (/ = 1) we have: 

^(1) 



D/2 


COS 6, 


r' 


<R' 


R' 


COS 6, 


r' 


>R', 



Vl'' = eZ{l-va,)x ^p1 (30) 



where S is the angle between the vectors r' and R' and V2 = —Vi (a^ — * —ctz)- 
Then the dipole term is given by 

r' 
r' < R' 

V^^^ = -2eZva,cos5x\ ^f (31) 



and in the dipole approximation 

1 r r' 

V = 1/(°) + V^^^ = 2eZ 



X { ^ — va^ cos 6x1 -^/ 



r' < R' 



^, r'>R' 



^/2 



(32) 



2.3 C. Exact multipole decomposition. The case 6 = 

For the calculation of electron-positron pair production in heavy ion reactions the 
situation as 6 ^ is of particular interest because it is in this lowest impact regime 
that the failure of perturbation theory is greatest. Since the scale at which nuclear 
effects contribute for two interacting large nuclei (6 < 15 fm) is much smaller than the 
atomic scale for electrons and positrons {h/rrieC = 386 fm), the 6 = solution should 
be of interest as closely approximating the situation at the smallest impact parameter 
for which nuclear effects are not applicable (or for which we omit all nuclear reactions). 
Sample calculations show 6 = and 6 = 15 fm results to be < 1% apart [15, 17]. 

Let us consider the 6 = case without any approximation in I/7. Here we assume 
t > 0. To obtain appropriate expressions for t < one can make use of the symmetry 
relation Mj^lr, —t) = (— l)'M°(r, t). By definition we have [15, 17] in the case of Vi : 



By symmetry only m = is no n- vanishing and we immediately integrate over (f){p = 
r sin^, z = rcosO) : 

j-i \r'^ sm t^/7"^ + (r cos d — vtY\^''^ 
Let us define x = cos 9 and note that 1/7^ = 1 — w^ to rewrite 

r^d), ,, _ \/2t(2/ + 1) /-i P,(x)dx 



Mr'ir,t)^^-^^j ., "'^^^ , (35, 

K + ^-l + x2^/2 



-^2 J.2 yrj.. 



This expression is symmetric in 1/v and t/r and the integral over x can be carried 
out in closed form for each value of /. The authors of Refs. [15, 17] recently found the 
following simple exact result for the multipole moments of the 6 = case: 

.1. „ A. 



n^ 2Jtt(21 + 1) Pii-)Qii-), r<vt 



rv 



Qi{-)Pi{-), r>vt. 
V r 



Then 



^ , 2.71(21 + 1) [ PiHQii-), r < vt 

I ^^ Qii-)PiC-), r>vt. 

y V r 

We obtain for V2 the expression 



V2 = eZ{l + va,)J2i-l)%%9)^L^_ ^x v r (3 



^ , , 2J7T{21 + 1) \ Pi{-)Qi{-), r<vt 

)V(-lW(^)-^^- ^x<^ n r 

^ - QKVK-), r>vt. 

K V r 



1. Monopole approximation (/ = 0) 
We obtain: 

Q«(-), r>vt. 

By taking the limit v -^ 0,lim -Qo(~) = l? so that in the non-relativistic limit we have 
Vnon-rei = 2eZ/r, as it should be. 

2. Dipole approximation {1 = 1) 
In this case we have: 

vfUeZ,l-.„„.)yOT2^xH*<^ ''^■'"' (40) 



rv 



-gi(-), r>vt 
r V 



and Vn 



(1) 



Vl'\a. 



-az). Then 



yd) = vl'^ + V, 



^(1) 



6eZa, 



cos^x <^ ^ 5; 



1^ /^ 
-Qi(-), r<vt 



Qii-), r>vt. 
r V 



Therefore in the dipole approximation 



V = \/(°) + \/(^) 



2eZ 




Sa^ cos 9 X < y V 



1^ ,t^ 
-Qii-), r<vt 



-Qii-), r>vt 
r V 



(41) 



(42) 



2.4 D. The case of large values of 7 and 6 = 

For 6 = the expression of multipole moments takes on the simple form good up to 
order In 7/7^ [15, 17]. For Vi : 



,-,, 2J7r(2/ + l) 



Qi 



r <t 



PK-)(ln27-^-), r>t. 



(43) 



n=l 



n 



Here t > 0. To obtain appropriate expressions for t < one can make use again of 
the symmetry relation. Analogously, we have for V2 '■ Mi {r,t) = (— 1)'M™ {r,t). 
In this asymptotic form, the split between the two regions r > vt,r < vt has been 
written as r > t,r < t, dropping the 1/7^ differences (we remind that w — > 1). 

1. Monopole approximation (/ = 0) 

We have: 

t. 



VI 



(0) 



eZ(l — a. 



In 27, r > t 



and W°^ 



V^/^«. 



—ay). Then 



(44) 



^(0) 



y(0) ^ yW + p. 



2. Dipole approximation (/ = 1) 
We have: 



2eZ 



X 



Qo 



r <t 



In 27, r > t. 



(45) 



^(1) 



\//^^ = eZ(l - a.)FA^) 



2V3^ 



X 



Qi(-), r<t 
r 

-(ln27-l), r>t 



(46) 
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and Vi^'^ = -Vl^\a, -^ -a,). Then 




Vii) = _^^^,osex{ , ^''r- ■ -" (47) 



Therefore in the dipole approximation 

2eZ 



V = \/(°) + \/(^) 



x^^o(^) -3«.cos.x| , ^^^^ '^' 



r 
In 27 I -(ln27-l), r>t 



■ (48) 



r 

In the monopole approximation, where only the term with / = in the muhipole 
expansion of the two-center potential is taken into account, a charged spherical shell 
of radius vt (case C) or t (case D) simulates the Coulomb potential of the two nuclei. 

3 Gauge transformation of the interaction 

It has previously been shown that in the large 7 limit a gauge transformation on the 
considered potential can remove both the large positive and negative time contribu- 
tions of the potential as well as the 7 dependence [15, 17]. If one makes the gauge 
transformation on the wave function [15, 17] 

i) = e-^'^^'^'^V, (49) 

where 



eZ 

V 

the interaction V{p,z,t) is gauge transformed to 



X{r, t) = — \n[-i{z - vt) + y^62 + ^2(^ _ ^^)2]^ (59) 



y (p, ., t) = , ^Z(l-va^) _ eZ(l-(l/v)a.) 

^[{h-p)hf+{z-vtY y^62/^2 + (^ _ ^^)2 

This transformation removes both the 7 dependence and the In 6 dependence of the 
interaction for the region r < V&^ + 1"^- The large 7 expression for the tti = gauge 
transformed interaction then becomes [15, 17] 

M^{r,t) = ^^ -Pi{-)x{ r'-e /-^-^ (52) 

r r -In — - — , r>Vb^ + t^. 

K b^ 

Note that this transformation affects the 771 = part of the interaction only. An 
additional computational advantage of the previous form Eq. (52) is that the large 

11 



positive and negative time contributions inherent in the term Qi{t/r) of Eq. (11) 
have been removed, as was shown by Baltz, Rhoades-Brown and Weneser [15]. It is 
very important that there is no 7 dependence in either m = or m ^ parts of the 
interaction. 

In the untransformed interaction, Eq. (11), the large negative and positive time 
interaction is dominated by the m = interaction for t > r, 



V,\r,t)=Y,\e)2V^^^?^Qiit/r) 



(53) 



(recall that there is no t > r interaction for m ^ 0). There is a similar Qi{t/r) 
dependence in the exact b = expression [15]. For the monopole interaction 



n. X 1 , 1 +r/t 
\/0 r,t = — In ^ ' 



I /j"^ /y^ 



(54) 



By analyzing Qi{t/r) for / 7^ 0, one finds that other 771 = interactions, V^°(r, t) have 
terms of longest range going as -7-7. For r < t we must consider Qi(t/r) which has a 
piece equal to — ln(t — r). Likewise, for t < r < ^/b"^ + t^, we must consider [15] 



In 



^.2 _-f-2 



62 



r + 1 
In —T7^ h ln(r — t). 



52 



1. Monopole approximatio7i 

In the case / = 0, m = 0, we have 



(55) 



Ko) 



1 



eZ(l-a^)— X 
Ir 



and W^ 



Vi («2 -^ —ctz)- Then 



0, r < Jb^/4 + t 



In 



62/4 



r > J62/4 + ^2 



(56) 



yiO) = vl'^ + V, 



.(0) 



2eZ 



2. Dipole approxi7natio7i (I = 1) 
The dipole term is given by 



0, r < J62/4 + t2 



In- 



r2-t2 
62/4 



r > JbyA + t^. 



(57) 



v;(i) = y/^^ + y, 



^(1) 



3eZa^ 



[sin ^ cos sin 6t x < 



' 62/4 ■ 
, 62/4 . 
'r2-t2- 



r <t 



t<r< Jh^/A + 12 



r > ./62/4 + t2 
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0, r < ^62/4 + ^2 
+ cos^x<( r2-t2 1. (58) 

In the dipole approximation we have to add the terms V^^^ and Vj-^^ given by Eqs. 
(57) and (58): 

K, = v;(°) + v;«. (59) 
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